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Abstract: We present the exact next-to- leading order (NLO) QCD corrections to di- 
jet production at the LHC via quark contact interactions, with different color and chiral 
structures induced from new physics. Following the recent analysis of quark compositeness 
search at the LHC, we find that the NLO QCD corrections can lower the dijet cross sections 
by several tens percent, depending on the theory parameters and the selected kinematic re- 
gions, and reduce the dependence of the cross sections on factorization and renormalization 
scales. We also calculate the renormalization group (RG) improved NLO cross sections by 
summing over the large logarithms from the RG running of Wilson coefficients. Moreover, 
we investigate the NLO QCD effects on various experimental observables and exclusion 
limits of quark compositeness scale. 
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1 Introduction 

Jet production at hadron colliders provides an excellent opportunity to test perturbative 
QCD (PQCD) and to search for possible new physics (NP) beyond the Standard Model 
(SM) over a wide range of energy scales. Invariant mass distributions of the dijets [1], the 
dijet angular distributions [2-4], and other jet observables at the LHC [5] have already 
extended current searches for quark compositeness, excited quarks, and other new particle 
resonances toward the highest energies attainable. Among all these measurements, the 
dijet angular distribution shows a great sensitivity to possible quark contact interactions 
induced by new physics models. In the SM, Quantum Chromodynamics (QCD) predicts 
the jets in dijet events are preferably produced in large rapidity region, via small angle 
scattering in t-channel processes. On the contrary, the dijet angular distribution induced 
by quark contact interactions is expected to be much more isotropic, and thus the dijet 
angular distribution at the LHC could be largely modified. 

The measurement of quark contact interactions has been used to set limits on the quark 
composite models which have been studied extensively in the literature [6, 7]. It is assumed 
that quarks are composed of more fundamental particles with new strong interactions at a 
compositeness scale A, much greater than the quark mass scales. At the energy well below 
A, quark contact interactions are induced by the underlying strong dynamics, and yield 
observable signals at hadron colliders. The newest bounds of A at the 95% confidence level 
(C.L.) from the CMS collaboration are around lOTeV [4] based on 2.2 fb~^ collected data. 
Previous limits from the Tevatron and LHC can be found in Refs. [2, 3, 8, 9]. With more 
integrated luminosity collected, these limits will be further improved. 
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In our previous work [10], we carried out the next-to-leading order (NLO) QCD cor- 
rection for dijet production at the LHC induced by the quark contact interactions that 
are the products of left-handed electroweak isoscalar quark currents. We compared it with 
the leading order (LO) results used by the CMS Collaboration and also the "scaled NLO 
results" used by the ATLAS Collaboration, which assumes the NLO correction (in terms 
of K-factors, defined as the ratio of NLO cross sections to LO ones) to the dijet production 
from contact interactions to be exactly the same as that from the SM QCD interactions. 
And we derived, based on our exact NLO results, the corrected limits of compositeness 
scale for the CMS and ATLAS measurements with 3pb~^ data [9]. In this paper, we ex- 
tend our previous work to include more quark contact interaction operators in the NLO 
QCD calculations, which allows the mixing of operators with different chiral structures 
at the NLO level, and show more details of the calculations. The effect of our results to 
measurements of quark compositeness at the LHC is also investigated. In the appendix, 
we discuss some details of our numerical code developed for the calculations presented in 
this paper. 

2 Theoretical setup 

We consider a subset of quark contact interactions that are the products of electroweak 
isoscalar quark currents which are assumed to be flavor-symmetric to avoid large flavor- 
changing neutral-current interactions [7]. The effective Lagrangian can be written as 



i=l 

where A is the new physics scale, Cj are Wilson coefficients. And the operators Oj in chiral 
basis are given by 




(2.1) 





(2.2) 
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in which c, d are generation indices and i, j, k, I, a are color indices, and are the 
GeU-Mann matrices with the normahzation Tr(T"T'') = 5°'^ /2. Beside of the quark com- 
positeness, the above interactions can also arise from various kinds of new physics models, 
induced by the exchange of new heavy resonances, such as Z' models [11] and extra di- 
mensions models [12]. Thus our analyses here are rather model independent and A can be 
identified as the effective new physics (NP) scale. 

The above six operators have been extensively studied in weak decays of mesons [13]. 
They mix with each other through QCD loop diagrams, which requires a renormalization 



matrix of the operators to cancel all the ultraviolet divergences, defined by O. 
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(1 + 



6Z)ijOj. After calculating one- loop diagrams, to be shown latter, with the dimensional 
regularization scheme in n = 4 — 2e dimensions, we obtain the matrix at the NLO as follow 
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(2.3) 



N 



3 and Cp = 4/3 for QCD, and nj = 5 is the activate quark 



numbers in the loop. (Here, we do not include the top quark contribution in the loops). 



3 Analytical results 
3.1 NLO corrections 

At LO, there are several subprocesses which contribute to the dijet production at hadron 
colliders induced by the operators under consideration. They are 

qq'{q) qq'{q), qq' qq', qq qq{q'q'), (3.1) 

where q, q' could be all the light quarks except the top quark. The NP contributions 
included in our calculation consist of two parts, the NP squared terms and the interference 
terms between the NP and the SM QCD interactions, which have different behavior with 
the increase of dijet invariant mass. We carried out the NLO calculations in the Feynman-'t 
Hooft gauge with dimensional regularization (DR) scheme (with naive 75 prescription) [13] 
in n = 4 — 2e dimensions to regularize all the divergences. Below, we only show the 
analytical results for the subprocess q{pi)q'{p2) <l{P'i)<l' {Pi)-, since the similar results for 
other subprocesses can be obtained by crossing symmetry. 

First, we define the following abbreviations for color structures and matrix elements. 

Ml = ULip3hiJ.ULiPl)uLiP4h^ULip2), 
M2 = UR{p3)^f,Uji{pi)uR{p4)-f^'uR{p2), 
M3 = {uL{p3)ltJiUL{Pl)uR{p4)l''UR{p2) 

+UR{p3)-1^UR{pi)UL{pi)-i^UL{p2))/'2, 
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Figure 1. Loop diagrams for botlr tire SM QCD and NP contributions to qq' qq' at the NLO. 
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(3.2) 



where ii_4 are the color indices of the external quarks. The LO scattering amplitudes 
induced by the NP and the SM QCD interactions can be separately written as 
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(3.3) 



where s, t, u are the Mandelstam variables, and we divide both the NP and SM QCD 
amplitudes into 3 groups for convenience. After adding the 1-loop amplitudes, as shown in 
Fig. 1, and the counterterms from renormalization, we obtain the ultraviolet finite virtual 
amplitudes as follows. 
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(3.5) 



The superscript uv in MJ^p is to indicate that they are ultraviolet finite. We have checked 
that the virtual correction for the SM QCD contributions given in Eq. (3.4) agrees with 
the ones shown in Ref. [14]. The infrared divergences in virtual corrections should cancel 
with those in real corrections. As for the real corrections, we apply both the phase space 
slicing based two cutoff method [15] and the subtraction based dipole method [16] in our 
calculations for a cross-check. The real emission diagrams include those shown in Fig. 2 
and all their crossing diagrams. 

For the new physics contributions, beside from the loop diagrams in Fig. 1, there are 
some additional penguin-like loop diagrams as shown in Fig. 3. They generate infrared finite 
terms with mixing among operators with different chiral structures. After renormalization, 
we obtain both infrared and ultraviolet finite amplitudes 



3(ci-^)(ln( — ) 
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Figure 2. Real emission diagrams for both the SM QCD and NP contributions to qq' qq' at the 
NLO. 



Figure 3. Additional penguin-like loop diagrams contributed to qq' — > qq' at the NLO. 
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which will contribute to the NLO results through interference with iM^^^p and i-M^J^j in 
Eq. (3.3). 

3.2 Renormalization group running of Wilson coefficients 

If the new physics scale A at which the Wilson coefficients are derived is much higher than 
the physics scale considered at colliders, there will be large logarithm terms associated 
with these two scales in fixed order calculations. We can improve the convergence of the 
perturbative calculation by summing the logarithm contributions using the renormalization 
group (RG) evolution of the Wilson coefficients. The RG equation is given by 



d In/XR 
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with the one-loop anomalous dimension matrix derived from Eq. (2.3), 
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The RG improved NLO cross sections are defined as 

O'NLCRG = O-NLO + (<7i?G " (^Rg\nLo) 







3N 

3Cf 
N 



(3.8) 



(3.9) 



where aj^c is the LO cross section calculated with RG improved Wilson coefficients and 
(^rg\nlo is the expansion of aRc up to NLO in QCD, which has already been included in 
o'NLO- If only ci or C2 is non-zero at the NP scale A, then the numerical solution of the 
RG equation that sums the LO logarithms are 



c^{^J■R) = ci{K)Xijr"<\ or Ci(/iij) = C2{h)Yijr"<\ 



(3.10) 



where r = as{fJ.R)/as{A). The eigenvalues of the one- loop anomalous dimension matrix 7, 
after being divided by the overall factor — ^ and /3o, are 



7. = (-0.942, 0.818, -0.627, 0.265, 0.183, -0.160), 



(3.11) 



where /3o = (llA^ — 2nj)/3. Furthermore, 
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(3.12) 



In Fig. 4 we compare fixed order (NLO) and RG running of the Wilson coefficients ci, C2 
and C4 for A = 30 TeV assuming only c\ (left panel) or C2 (right panel) is non-zero at the 
NP scale A, with the renormalization scale ranges between 1 TeV and 30 TeV. The running 
of the Wilson coefficients at the NLO can be obtained from the expansion of RG running 
in Eq. (3.10). 
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HR(TeV) ^^(TeV) 

Figure 4. NLO and RG improved running of the Wilson coefficients. Curves from top to bottom 
correspond to ci (black), C4 (blue) and C2 (red) for the left panel, and C2 (red), ci (black) and C4 
(blue) for the right panel. 

4 Applications to quark compositeness search at the LHC 

In this section we apply our NLO QCD results to the quark compositeness search at the 
LHC {^/s = 7TeV) through dijet angular distribution measurement. For the numerical 
results here we assume only ci(A) and C2(A) to be non-zero, and parameterize them as 
Ci(2)(^) = 47rAi(2)- Following conventions used in experimental analysis, for color-singlet 
case we have |Ai(2)| = 1(0), while for color-octet case we have |Ai(2)| = 0(1). We use the 
anti-fcj jet algorithm [17] with energy recombination scheme [18] and the distance parameter 
D = 0.5. To be considered as one of the two leading jets, a jet is required to satisfy the 
rapidity cut \y\ < 3. Moreover, we apply additional constraints on jet rapidity as below 

\yb\ = \yi + y2|/2 < i.n, x = exp {\yi - 2/2I) < I6, (4.1) 

to be consistent with the CMS measurement [2]. x is chosen as the dijet angular observable 
since after the Jacobian transformation the t-channel dominant SM QCD dijet production 
is almost flat on x- For massless jets, the invariant mass of dijet can be expressed as 

rrijj = y/PTiPT2Vx + '^/x - 2cos(A(/>), (4.2) 

where A(f) is the azimuthal angle between the two jets. Considering the current experi- 
mental limits on the compositeness scale as well as the large SM QCD dijet background, 
NP contributions can only provide observable effects on the angular distribution in a very 
high dijet invariant mass region. Thus, below we will consider an invariant mass region 
between [2TeV and 3TeV], and the part of theory parameter space with A greater than 
5TeV, for simplicity. In our numerical calculations, we use CTEQ6.6 parton distribution 
functions [19] and the corresponding running QCD coupling constant. Renormalization 
and factorization scales are set to be the average transverse momentum of the two leading 
jets, unless otherwise specified. 
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4.1 LO results and analysis 

The LO total cross sections for dijet production induced by contact interactions consists 
of interference contributions with SM QCD amphtudes (denoted as ctint), as well as the 
NP squared contributions (denoted as cfsq)- Assuming that only ci and C2 are non-zero, 
the additional LO contribution 

<7LO = <yiNT + cySQ (4.3) 

= (Ai6l,i + A26l,2)/A^ + (A?6l,ii + A^6l,22 + AlA26L,l2)/A^ 

where hi^^n^ij-^ are independent of NP scale A and c, = 47rAj. Instead of calculating the 
differential cross sections with respect to X) we choose two representative bins to investigate 
the influence of NP contributions to the dijet angular distribution, which are x = [1)6] (bin 
1) and X = [6, H] (bin 2). The numerical results of bL,j(ij) are shown in Table 1. It can be 
seen that absolute values of h for bin 1 are much larger than those for bin 2, which means 
the NP contributions could change the shape of x distribution significantly since the SM 
QCD contributions are almost flat in x distribution. It is also indicated in Table 1 that, 
even for a scale of A as large as 5TeV, the NP squared terms can still be comparable to 
or even larger than the interference terms, especially for bin 1 of color-singlet case. Thus, 
the cross sections are not monotonous decreasing as A increases, and zero point occurs for 
certain A value with destructive interference. For simplicity, we will focus on the parameter 
region where the NP squared terms are relatively small, i.e., (JsQ/f^iNT < 1. Fig. 5 shows 
the absolute ratios of the NP contributions from the squared terms to the ones from the 
interference as functions of A, where for the color-singlet (octet) case we set |Ai| = 1 (0) 
and IA2I = 0(1) following the standard convention as mentioned before. Using the above 
condition, constraints on A would be > 8TeV for the color-singlet case and > 4TeV for 
the color-octet case, which will be applied hereafter to our analysis. 



[fb-(5TeV)2(4)] 




bL,2 


feL,ll 


&L,22 


&L,12 


bin 1 


-258 


-179 


614 


93.4 


259 


bin 2 


-99.1 


-70.4 


113 


17.2 


46.8 



Table 1. LO coefficients of NP contributions. 



4.2 NLO results and analysis 

At NLO, the total cross sections for dijet production induced by the contact interactions 
can be expressed as 

c^NLO = (Ai(6n,i + air) + A2(6n,2 + a2r))/A^ + (A?(&n,ii + onr) (4.4) 
+A2(^'N,22 + a22r) + AiA2(6n,12 + ai2r))/A^, 

where &N,j(ij) and a^^ij) are independent of the NP scale A, and r = ln(A/po). We choose 
the reference scale 

Po = {mjj)/V{x) + l/ix) + 2, (4.5) 
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Figure 5. Ratios of NP squared contributions {(Tsq) to SM and NP interference contributions 
(ctint) at tlie LO. 

where (x) denotes taking the average value of x in the given bin of kinematic variables. 
Hence, po is 1.04 (0.77) TeV for bin 1 (2). 

The contributions proportional to r represent running effects of the Wilson coefficients 
at NLO. Choosing same bins as in the LO analysis, we get the results for iN,i(ij) &iid ai[ij) 
listed in Table 2. Based on results from Tables 1 and 2, we plot the NLO K- factors as 
functions of the compositeness scale in Fig. 6 for both cases with destructive and construc- 
tive interferences. The shadow regions indicate parameter spaces where contributions from 
the NP squared terms are larger than the ones from the interference terms. K-factors are 
unstable there for destructive interference case due to large cancelations of cross sections 
between NP squared contributions and interference ones. Beyond the shadow regions, the 
NLO QCD corrections reduce the absolute values of the NP contributions to the cross sec- 
tions significantly for all the parameter cases. These are mainly due to the large constant 
terms and also the large logarithms of A from the virtual corrections, cf. Eqs. (3.4) and 
(3.5) in Sec. 3.1, especially in the large x region. From Fig. 6 we can see that the K-factors 
deviate significantly from 1 especially for bin 2. Thus one may doubt the reliability of the 
NLO results, i.e., the convergency of the perturbation series. Indeed, the small K-factors 
are mainly due to the LO results used here. In order to compare with the LO theoretical 
results used by the experimentalist, by default we use fixed LO Wilson coefficients in the 
LO cross sections here and below. If we use the NLO running Wilson coefficients in the LO 
calculations instead, same as in the NLO calculations, then the K-factors will increase due 
to the suppression of the LO cross sections. Also the K-factors depend on the QCD scale 
choices, especially for large x region. If we set the central scale to PT,max exp(0.15|yi — ^2!) 
as in the ATLAS dijet study [20], where PT,max is the transverse momentum of the hardest 
jet, then the K-factors will further increase. In Fig. 7 we show the K-factors calculated 
using this alternative definition of the LO cross sections in the perturbation series and 
choice of the scale. We can see that the NLO results are well behaved with reasonable 
K-factors that are relatively larger and stable as compared to the ones in Fig. 6. 
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Table 2. NLO coefficients of NP contributions. 




A (TeV) A (TeV) 

Figure 6. NLO K-factors, as functions of A, using the default definition of the LO cross sections 
and choice of scale. 




A (TeV) A (TeV) 

Figure 7. NLO K-factors, as functions of A, using the alternative definition of the LO cross sections 
and choice of scale. 

In Fig. 8 we show NP contributions to the total cross sections of bin 1 with scale un- 
certainties to further investigate the improvement on scale dependence of the NLO results. 
The uncertainties are calculated by varying the factorization and renormalization scales, 
separately, for fj, = (px) /2, fi = (pt), and = 2{pt), where (pt) is the average pT of the 
two leading jets. Generally, we can see a reduction of the scale uncertainties for the NLO 
results. 

4.3 RG improved NLO results 

As already mentioned in the analytical results, when the compositeness scale is much 
higher than the typical energy scale of experiments, we can sum over large logarithms 
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Figure 8. LO and NLO total cross sections, including scale uncertainties, as functions of A. 



of A, induced from higher order corrections, to improve the convergency of perturbative 
expansion. This leads to the RG improved NLO cross sections ctrg.nlo- It can deviate 
largely from the NLO cross sections in the large x region where the jet pT is relatively 
small so that the running effects of Wilson coefficients are important. In Fig. 9, We show 
the ratios o'rg,nlo/'7lo and (Tr,g,nlo/c"nlo as a function of A. In general, the higher order 
corrections can increase the NLO total cross sections for both color singlet and octet cases, 
and the amount depends on the kinematic region considered. For example, for the bin 1 of 
color singlet case, the increase in ctrq^nlo is about 3% and 10% of ctnlo for A = 8 TeV and 
30 TeV, respectively. Moreover, the resumed contributions stabilize the K-factors of the 
total cross sections as compared to the NLO results for high A values, which are around 
0.7 and 0.5 for the bin 1 and the bin 2 of the color singlet case, respectively. 





Figure 9. K-factors of the RGE improved NLO cross sections, as functions of A. 



4.4 Exclusion limits of quark compositeness scale at the LHC 

In order to directly compare our predictions to the experimental measurements, we need 
to combine them with the pure SM QCD contribution calculated at the NLO. We use 
a modified version of the EKS code [21] to calculate the SM QCD dijet production at 
the NLO. We only consider the case of color-singlet with constructive interference in this 
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Figure 10. Dijet invariant mass distributions from the pure SM QCD contribution at NLO and 
the NP induced ones at both LO and NLO. 



section, and we do not include the RG improved corrections here since they are small for 
the considered A values. In Fig. 10, we show comparisons of differential cross section from 
NP and pure SM contributions with invariant mass ranges from 2 to 3 TeV for the color- 
singlet case with constructive interference and A = 8 TeV. As already mentioned before, 
the pure SM contributions in the two bins are almost the same, while the NP contributions 
are greatly different. To derive the expected exclusion limits of the compositeness scale, we 
further divide the invariant mass region [2 TeV, 3 TeV] into 10 mass bins with equal width, 
and define the measure in each mass bin, F^{Mjj) = (7uni{Mjj)/aun2{Mjj). In Fig. 11, 
we plot theoretical predictions for F^{Mjj) in different mass bins, where the (gray) solid 
vertical line represents the pseudo data expected by pure SM QCD contributions. The 
errors of the pseudo data include both the estimated statistical and systematical errors. 
The former is estimated by assuming Gaussian statistics with one standard deviation. 
Here, we assume a large data set with C = 5 fb~^ to calculate the statistical errors. We 
note that most of the experimental systematic uncertainties cancel in the ratio F^{Mjj). 
For simplicity, we estimate an overall experimental systematic uncertainty of 3%, arising 
from the jet energy calibration and jet pT resolution for all the mass bins [2] . The (colored) 
dashed and solid horizonal lines represent the SM QCD contributions plus NP contributions 
at LO and NLO, respectively. Fig. 11 shows that the NLO QCD corrections reduce the 
NP contributions to F^{Mjj). 

On the other hand, the theoretical predictions also have uncertainties due to parton 
distribution functions, non-perturbative corrections, and most importantly the unknown 
higher order QCD corrections. The first two are found to be small according to the analysis 
in Ref [2]. The conventional way to estimate the uncertainty from unknown higher order 
corrections is to examine its scale variations, i.e, to calculate the spread of the cross section 
over a set of scale choices. In Fig. 12(a), we show scale variations of the cross sections 
in different bins for SM plus NLO NP contributions, the spread are calculated by varying 
the factorization and renormalization scales independently for fi = {pt)/2, ^ = (pt), and 
fi = 2{pt)- For F^{Mjj) = cruniiMjj) / aun2{Mjj) , we can independently vary the scales in 
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Figure 11. F^{Mjj) for pseudo data estimated by pure SM QCD contribution, and for SM plus 
NP contributions at the LO and NLO, respectively. 

both the numerator and denominator factors. Since there may be some correlations between 
these two parts in the missing higher order corrections, it will lead to an overestimation of 
the uncertainties. Here, we vary the scale simultaneously for cross sections in bin 1 and bin 
2, and take half of the total scale variation of F^{Mjj) as the estimated theoretical error. 
The corresponding theoretical errors of Fy.{Mjj) for SM plus NLO NP contributions are 
shown in Fig. 12(b), which are around 4-7% in different invariant mass bins. 
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Figure 12. (a) Total cross sections in bin 1 and bin 2 with scale variations; (b) predictions of 
F^{Mjj) with theoretical errors included. Both of them are for SM plus NLO NP contributions. 



We perform a simple log-likelihood t^st for the hypothesis of NP with 
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1=1.10 



(4.6) 
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where F^^{i) represents the pure SM contribution of F-^{Mjj) in the ith mass bin, which we 
assume to be the expected data, and F^^^~^^^{i) is the theory prediction given by SM plus 
NP contributions. Ai.xp,th represents the corresponding experimental errors and theoretical 
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errors of F^^^^'^^{i), and we do not consider possible correlations of errors in different mass 
bins. The x^/^d.o.f with N^.o.f = 10 are shown in Fig. 13 as functions of the compositeness 
scale A for 3 cases, i.e, SM plus LO or NLO NP contributions without including the 
theoretical uncertainty, and SM plus NLO NP contributions with theoretical uncertainty 
included. Exclusion limits (95% C.L.) of A can be read directly from Fig. 13 as intersections 
of the curves with the horizontal line. We can see that the theoretical uncertainty has a 
large effect on the exclusion limit since they are comparable with the experimental errors. 
With more data collected at the LHC, the statistical error will further decrease, then 
theoretical uncertainty will play a much more important role in the measurement. 
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Figure 13. /^d.o.f as functions of the compositeness scale A. Intersections of the curves with 
the horizontal line show the exclusion limits at the 95% C.L.. 



5 Conclusion 

In conclusion, we have calculated the exact NLO QCD corrections to dijet production at the 
LHC, induced by quark contact interactions with different color and chiral structures from 
new physics. By applying our results to quark compositeness search at the LHC, we show 
that the NLO QCD corrections can lower the NP cross sections by several tens percent, 
depending on the choice of the theory parameters and the kinematic regions considered. 
Moreover, the NLO QCD corrections reduce the dependence of the cross sections on factor- 
ization and renormalization scales. We also calculate the renormalization group improved 
NLO cross sections by summing over the large logarithms induced from the running of 
Wilson coefficients, which are found to stabilize the K-factors at large quark compositeness 
scales. We further investigate the NLO QCD effects on the corresponding experimental 
observables and study the exclusion limits of the quark compositeness scale. 

Appendix 

In this appendix we briefly introduce our numerical code CIDIJET2.3, developed for the 
NLO QCD calculation of the dijet production induced by quark contact interaction. 
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Basically, this code is calculating the double differential cross sections in terms of the 
dijet invariant mass Mjj and the angle parameter x- There are two calculation modes in 
this code. The first one is to numerically evaluate the fixed order double differential cross 
sections, in a kinematic bin specified by its range of Mjj and for various values of A and 
Cj (or Aj). For Cj = 47rAi, the fixed order (LO or NLO) cross section can be written as 

6 6 
i=l i=l 

+ ^ {\i\i+i{hii+i + au+ir)) / 

1=1,3,5 

+ i>^MbiA + ai4r))/A\ (5.1) 

j=l, 2,5,6 

with r = ln(A/po), and po is defined as in Eq. (4.5). The code can calculate all the above 
coefficients (a's and 6's) directly, thus users do not need to repeat the calculations for 
different A and Aj values. Since QCD interaction preserves parity conservation, we should 

have 6(a)i,2 = &(a)5,6, b(a)ii,22,i2 = ^(0)55,66,56, and 6(0)14,24 = 6(0)54,64- The code has 
two dynamic QCD scale choices, which is average pT of two leading jets used by the CMS 
collaboration and PT,max exp(0.15|2/i — 2/2!) used by the ATLAS collaboration, and provides 
inputs of pre-factors for both renormalization and factorization scale to study the scale 
variations. Another mode is for the calculation of higher order corrections in addition to 
the NLO cross sections, which arise from the renormalization group running of the Wilson 
coefficients, as discussed in section 3.2. They are only significant at very large A values and 
need to be recalculated for different inputs of A and Aj. The CIDIJET2.3 code is publicly 
available by request. 
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